Trapped fermionic clouds distorted from the trap shape due to many-body effects 
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We present a general approach for describing trapped Fermi gases, when the cloud shape is 
distorted with respect to the trap shape. Our approach provides a consistent way to explore physics 
beyond the local density approximation, if this is necessary due to the distortion. We illustrate this 
by analyzing in detail experimentally observed distortions in a trapped imbalanced Fermi mixture. 
In particular, we demonstrate in that case dramatic deviations from ellipsoidal cloud shapes arising 
from the competition between surface and bulk energies. 
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Introduction. - The cooling of trapped fermionic 
atoms to degeneracy has opened up a fascinating new 
arena for studying many-fermion physics In partic- 
ular, the opportunity to access the BEC-BCS crossover 
using Feshbach resonances has led to a large effort in the 
experimental study of fermionic pairing using trapped 
atomic gases [2j. Among other topics, pairing in polar- 
ized Fermi gases is now being intensely explored, experi- 
mentally in Refs. [a, \4 , la [ql and theoretically in, e.g., 
Refs. 0, i 0, [lliHllT Apart from possible new 

insights for solid-state and nuclear physics, the field of 
trapped fermionic atoms comes with its own set of new 
concerns, such as the role of the trapping potential. In 
this Letter, we therefore provide results on the relation- 
ship between the trap shape and the shape of the gas 
cloud. 

For a large class of trapped many-body systems, the 
local density approximation (LDA) gives an accurate de- 
scription of the spatial distributions of atomic and energy 
densities. However, motivated by recent experiments 
0, , we propose that a number of mechanisms may lead 
to the intriguing situation that the combination of many- 
body physics and the external trapping potential causes 
the shape of the gas cloud to be distorted compared to the 
shape determined by the trap alone. In such cases, the 
simple LDA procedure in terms of a local chemical po- 
tential is no longer appropriate, and we need to consider 
the angular structure of the local Fermi surface. We have 
therefore formulated an extension of the LDA that allows 
for a distinction between the different directions. We give 
below some examples where one observes or expects such 
spatial distortions. We note however that there might be 
many other situations, beyond what we anticipate here, 
where such distortions take place. Our formalism is thus 
of very general applicability and interest. 

As a first scenario for spatial distortions, consider a 
many-body system that is phase separated in an elon- 
gated trap due to the presence of a first-order transi- 
tion. The interface between the two phases then carries 
a surface tension. To minimize the surface energy, the 
interface shape will differ from the trap shape. In fact, 
such a system has recently been experimentally realized 



using a trapped polarized Fermi mixture of resonantly- 
interacting atoms 0, Q . At sufficiently low temperatures 
this mixture phase separates into an unpolarized super- 
fluid core and an outer shell of normal polarized fermions 
0, 0, H @, EU EH- When the trap is elongated enough 
and the particle number is small enough, the surface ten- 
sion causes significant deviation of the core aspect ratio 
from the trap aspect ratio [f| [l2 |. In the following we 
will analyze this situation in detail. 

Second, Pomeranchuk instabilities are shape deforma- 
tion instabilities of a Fermi surface [T3|, which have at- 
tracted intense interest in the strongly correlated elec- 
trons community [141 1 1 51 ] . A <i-wave Pomeranchuk tran- 
sition leads to a "nematic" Fermi liquid with an ellip- 
tical Fermi surface. Other than their possible realiza- 
tion in solids, it is hoped that such a nematic liquid may 
be realized in atomic dipolar Fermi gases since a long- 
range interaction can cause a Pomeranchuk instability if 
the short-range part is suppressed [l5j |. A Fermi surface 
deformation that is constant in space causes no distor- 
tion in coordinate space. However, in an external trap, 
the Fermi surface deformation would vary from point to 
point, leading to a distortion of the gas cloud in coordi- 
nate space as well. To calculate particle and energy den- 
sities in such a case, would again require an anisotropic 
generalization of the LDA, of the kind we provide here. 

Finally, one possible mechanism for pairing of polar- 
ized fermions is via opposite deformation of the two Fermi 
surfaces [3]. A prolate (oblate) shape for the major- 
ity (minority) Fermi surface leads to an equatorial re- 
gion in momentum space along which opposite-momenta 
fermions can pair. With the extended parameter space 
expected to be accessible with trapped atoms in the fu- 
ture, it is likely that this pairing mechanism is energeti- 
cally favorable in some parameter regime [r| . In the case 
of such pairing in a trap, we expect both Fermi surfaces 
to have ellipticities (of opposite sign) that vary with po- 
sition. As a result, the majority and minority clouds are 
expected to be distorted in opposite directions relative 
to the trap. Again, a treatment of the type presented in 
this Letter is needed to describe this situation. 

Anisotropic generalization of local density approxima- 
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Hon — The LDA encodes the inhomogeneity of the gas 
completely in terms of a local chemical potential at ev- 
ery point in space, //(r) = \x — V tr (r). Here \x is the true 
chemical potential of the gas and V tI (r) is the trapping 
potential. The shape of the gas cloud is determined by 
an equipotential surface of the trap. For harmonic traps, 
this surface is an ellipsoid. 

The key to treating distorted systems is to general- 
ize the definition of the local chemical potential /i(r) = 
A* - Vtr(r) in a way that discriminates between directions. 
For definiteness, we will consider axially symmetric traps, 
Vtr(r) = V tI (r,z), where r and z are the radial and ax- 
ial length variables of a cylindrical coordinate system, 
respectively. For each position in the trap, we need to 
determine the momentum-space location of the deformed 
Fermi surface, in both the radial and axial directions. We 
denote by R(z) and Z(r) the edge of a reference ideal 
Fermi cloud, that would be obtained by hypothetically 
turning off the interactions but retaining the distortion. 
Our prescriptions for the location of the Fermi surface at 
position r in the radial and axial directions are 

A*R.( r ) = V>B.(r,z) = Vtr(R(z),z) -Vtr(r,z) , (la) 
/i Z (r) = n z (r,z) = V tI (r,Z(r))-V tI (r,z) . (lb) 

We have written the Fermi surface locations in energy 
units, i.e., /xr and \xz are the free- particle energies corre- 
sponding to the Fermi momenta in the radial and axial 
directions. They are, of course, not components of the 
chemical potential. Eqs. fl} are obtained by semiclassical 
considerations similar to the WKB approximation, which 
is one way to derive the LDA. Here, we have applied the 
semiclassical approximation separately to the r and z di- 
rections (Fig. [TJi) as appropriate for an axially symmetric 
trap, where the single-particle Schrodinger equation sep- 
arates in these two directions. 

After obtaining the locations of the local Fermi surface 
in the radial and axial directions, it is for many purposes 
sufficient to use an ellipsoidal Fermi surface with /zr and 
/iz as the principal axes. At the center of the trap, the 
same semiclassical argument can be carried out in any 
direction and the local Fermi surface has exactly the same 
shape as the edge of the ideal gas reference cloud. The 
cloud shape is generally near ellipsoidal but may deviate 
from this shape either due to anharmonicity of the trap 
or due to many-body effects, as we will see. Although the 
shape of the Fermi surface away from the trap center is 
not fixed by similar arguments, it is physically reasonable 
to use the same shape as in the trap center. 

The local non-spherical Fermi surfaces obtained us- 
ing this prescription can now be used to calculate den- 
sities. For example, for an clliptically distorted normal 
gas in a harmonic trap, the density is given in terms 
of the volume of the local Fermi surface by n(r) = 
[(2m) 3 / 2 /67r 2 ] /zr(i")-\/ /iz(r). This density vanishes at 
the surface defined by R{z) or Z(r). For a superfluid, on 
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FIG. 1: (Color online.) (a) Calculation of the local 
anisotropic Fermi surface at position r (full dot). Points 1 
and 2 with coordinates (±R(z),z) are the classical turning 
points for the motion in the r-direction. Points 3 and 4, i.e., 
(r, ±Z(r)), are the classical turning points for motion in the z 
direction, (b) Fits to the axial density profiles for both species 
using our approach. The top and bottom panels show two 
experimental shots, at polarizations P ~ 0.5 and P ~ 0.65, 
respectively. 



the other hand, it may be appropriate to use the BCS 
equations generalized to anisotropic Fermi surfaces, i.e., 
the momentum integrand in the density profile 
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now has angular dependence. Here x = cos e k = 
fi 2 k 2 /2rn with m the atomic mass, A is the gap, and 
€f{x) is the deformed Fermi surface with principal axes 
/LtR and Similar modifications can be made for the 
BCS energy equation and the BCS gap equation. 

Application to polarized resonantly-interacting Fermi 
gases — We now demonstrate our procedure by elabo- 
rating on a particular case, i.e., we analyze the experi- 
ment reported in Ref. Q using the formalism presented 
above. Figure 2 of Ref. shows that the majority as- 
pect ratio does not change much with polarization. We 
therefore describe the majority species in state |f) with 
the usual LDA. On the other hand, the boundary of the 
minority atom cloud in state | J. ) is distorted due to sur- 
face tension. We thus treat the minority according to the 
approach described above. At each position within the 
minority cloud, we calculate the two parameters for the 
Fermi surface of the minority species, using as an input 
the boundary of the minority cloud, which coincides with 
the superfluid core if any intermediate partially-polarized 
normal shell can be neglected. 

For a choice of /ij, R^ = Ri(0), and Z^ = Zi(0), the 
particle densities n^(r) = n|(r) and the energy density 
e(r) at any point inside the core can be calculated from 
the local Fermi surfaces. The majority Fermi surface 
is given by /if(r) = //j — Vt r (r). The minority Fermi 
surface is deformed, with /XRj(r) and /zzj(r) calculated 
using Eqs. {1}. Since at resonance the gap is larger than 
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the maximum difference between the Fermi surfaces, the 
minority momentum distribution is not anisotropic even 
though fj, R ± ^ nz\ , unlike the case of Refs. [l6| • For the 
fully polarized normal gas outside the core, n^(r) and 
e(r) are calculated by the usual LDA. 

To calculate the densities and superfiuid gap in the su- 
pcrfiuid core, we could use generalized BCS equations as 
described above. However, with one additional approxi- 
mation, we can also make use of Monte Carlo results [13] 
for the resonant situation near which the experiments of 
Refs. are performed. The additional approximation 
is to use at every position an effective chemical potential 
/ i l( r ) — [M-R|( r )] 2//3 [Mzi( r )] 1 ^ 3 f° r the minority species 
as well. The superfiuid properties are then given by the 
universal unitarity results [17| in terms of the average 
chemical potential /i(r) = [/i|( r ) + Mt( r )]/2- 

In Fig.OJ), we show the resulting axial density profiles, 
and compare with experimental data for two different po- 
larizations P = (Nj - Ni)/(N^ + N±). The calculation 
requires three parameters (/if, R±, Z[) as input. Unlike 
the calculation in Ref . [1(3] , the required input quantities 
are not immediately fixed by the trap parameters and 
particle numbers (Af,Aj) alone. An additional energy 
calculation is required for the aspect ratio of the core, as 
explained later. The experimental majority axial density 
profiles have pronounced peaks at the center. This is re- 
produced well by the calculations shown in Fig. [T] using 
Monte Carlo unitarity parameters, but not if we use the 
BCS equations to describe the core. Note that these cen- 
tral peaks are not prominent in the first experiment by 
Partridge et al. [jjj], nor in the experiment by Shin et al. 
(H , suggesting that thermal fluctuations play a stronger 
role in those cases. In all calculations reported in this 
Letter, anharmonicities of the experimental trap shape 
are included unless explicitly stated otherwise. 

Energy minimization — At unitarity, the surface en- 
ergy density associated with the superfluid-normal inter- 
face is expected to be given by the chemical potentials 
through universal constants. Defining ft = (/x-j- + fJ>i)/2 
and h = (/if — fii)/2, we first note that a quantity with 
the dimension of a surface energy density constructed 
out of ft alone must be proportional to (m/h 2 )ft 2 . The 
dependence on h may be included as a function of the di- 
mensionless ratio h/ ft. The surface tension thus has the 
form a = rj s {mj 'h 2 )ft 2 f s (h/ ft) . Here r} s is a dimensionless 
number expected to be of order 1 and / s is a dimension- 
less function. Noting in addition that the position of the 
interface is fixed by the first-order transition determined 
by a universal value of h/ ft [l7l ]. we infer that f s (h/fi) 
is a constant along the surface we are interested in. We 
thus absorb it into the universal constant r/ s . 

We can use the observed core aspect ratios Z±/R± 
extracted from the axial density data to determine T} s , 
because the correct value of r/ s must produce a total 
energy minimum at the correct aspect ratio. Denot- 
ing by and e^(n) the particle and energy-density 
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FIG. 2: (Color online.) (a) Total energy curves for fixed 
(N h N L ) = (194,92.7) x 10 3 , shown for n s values of 0.30, 
0.59 and 0.90 from bottom to top. (b) Core aspect ratio 
(A.R.) as a function of polarization for iV T /10 3 = 200 - 80P 
and calculated with n s ~ 0.6. (c) Core A.R. versus particle 
number, for fixed polarization P = 0.49. The regions relevant 
to experiments 0, 0, [HUt are indicated by shaded ellipses. 



functions of a homogeneous two-component ideal Fermi 
gas, the bulk energy density in the superfiuid core is 
given by e(r) = (1 + f3)ej$ (n(r)) + V tl -(r)n(r), where 
ri(r) = % (/i(r)/(l + j3)) is the density profile in the core, 
[3 ~ —0.585 is a universal constant [lj], and ft(r) is cal- 
culated as described previously. The energy density in 
the outer shell is given simply by the single-component 
ideal Fermi gas expression. Integrating over the core and 
the outer shell gives the total bulk energy Et,. 

The total energy is Et, + i]sls, where I s is the surface 
integral over (mj h 2 ) [x 2 (r) . In Fig. [2k, this quantity is 
shown for several values of ?7 S , for fixed populations. The 
axial data in an experimental shot for these populations 
is fit well with Z^/R± ~ 21.5 and a minimum is located 
at this core aspect ratio for rj s ~ 0.59. Carrying out 
this procedure with several experimental axial density 
profiles, we obtain the estimate r; s = 0.60 ± 0.15. The 
error bar leads to uncertainties in the calculated aspect 
ratios that is about 2% for the smallest deformations and 
about 10% for the largest deformations. 

Variation of core deformation — Having determined 
the universal constant 7; s , we can compare our theory 
to the experimental core aspect ratio versus polarization 
curve in Figure 2 of Ref. [6j. The aspect ratio actually 
depends on the total particle number in addition to the 
polarization. In the experiment, there are shot-to-shot 
variations of N-\ and a systematic bias towards larger 
2Vf at smaller P. We model this bias as iV T = (200 — 
80P) x 10 3 . Most of the shots are scattered around this 
line. Our results are shown in Fig. [2Jd, together with the 
experimental data. 

In Fig. [2J;, we display the core deformation as a func- 
tion of total particle number. Surface effects are less 
important in the thermodynamic limit N — ► oo and the 
difference between core and trap aspect ratios vanishes 
as A -1 / 3 in this limit [12]. The lower curve is for a har- 
monic trap matching the trap of Ref. 0] at the center. 
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FIG. 3: (Color online.) Experimental axial density curves 
fitted with core shapes (a) 7 = 2 and (b) 7 = 8. Top panels 
show species axial densities, bottom panels show difference. 
The value 7 = 2 misses the "flat" features, while 7 = 8 over- 
emphasizes them, (c) Surface shapes for 7 = 2, 4, and 9, 
respectively, (d) Calculated value of 7 for various polariza- 
tions, along the line 7V T /10 3 = 200 - 80P. 



sity data because of the large fluctuations of (iVf, N\) 
from the line 7V T /10 3 = 200 - 80F. 

Conclusions — To summarize, we have presented an 
anisotropic extension of the local density approximation 
in order to deal with situations where the aspect ratio of 
trapped fermionic gases differs from the trap aspect ra- 
tio due to many-body effects. Our procedure allows for 
a consistent calculation of the density and energy pro- 
files. We have demonstrated this via a detailed analysis 
of experiments with phase-separated polarized fcrmions. 

In particular, we have uncovered and explained a strik- 
ing new effect concerning the shape of the superfluid core, 
which takes an unexpected non-ellipsoidal form. Our 
treatment dramatically highlights the utility of our for- 
mulation: we know of no other way of calculating the 
density distributions for an arbitrary cloud shape. 
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Hulet. 



The upper curve presents our expected result for the ex- 
periments of Shin et al. [5| . The predicted deformation of 
about 10% appears to be sufficiently large to have been 
observable in these experiments, which is consistent with 
the interpretation of Gubbels et al. [l8| that no phase 
separation has occurred in this case. 

Deviation from ellipsoidal shape — Examination of the 
high-polarization data of Figure 1 of Ref. [6j reveals a 
peculiar feature. The minority absorption images, other 
than having a lower aspect ratio than the trap, also have 
shapes that are not quite ellipsoidal and tend to being 
somewhat cylindrical. The axial density profiles reveal 
the same geometric effect in terms of flatter peaks than 
would be expected for an ellipsoidal shape. Our formal- 
ism is well-suited for studying this effect, since the core 
boundary simply translates into distortions of local Fermi 
surfaces that can be used to calculate consistent den- 
sity and energy profiles. We parameterize non-elliptical 
shapes by describing the core surface by the equation 
(r/i?) 7 + (z/Z) 7 = 1. Using a 7 equal to 2 gives an el- 
lipsoid and infinite 7 gives a cylinder. Comparing axial 
density profiles obtained with various values of 7, as in 
Figs. [3^ and Ob, shows that the experimental features, 
e.g. the flat tops, are best reproduced with 7 between 3 
and 6, with the best 7 tending to increase with increasing 
polarization, up to about P ~ 0.6. However, a precise 
determination of 7 is difficult from these fits alone. 

Our formalism also allows us to extract 7 from an en- 
ergy calculation. Minimizing the total (bulk plus surface) 
energy, under the constraint of constant particle num- 
bers, produces the observed deviation from ellipsoidal 
shape, with 7 > 2 and 7 increasing with P for moder- 
ate polarizations (Fig. [5Ji). The energy calculation also 
predicts a decrease of 7 at larger P. This prediction is 
difficult to verify from the currently available axial den- 
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